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Abstract 

In this paper, a new method based on signal flow graph technique and Mason's gain formula are applied to model a dc-dc Zeta 
converter. The model is the transfer functions from control to output and from input to output. Zeta converter is a rather 
complex system, which has four passive elements and may perform in either buck or boost operating modes. Compared with 
other converters in the same group, such as Cuk and SEPIC converters, the dynamic modelling and analysis of the Zeta 
converter have been not reported in the literature. For modelling target and to get rid of complex equations in modelling, a new 
method based on signal flow graph approach is applied to manipulate the equations. After modelling and doing dynamic 
stability analysis in continuous conduction mode, and prove the necessity of a controller design for it. A model based PID 
controller is designed for the converter and the simulation results are presented to verify the correctness of the achieved model 
and model- based controller and the signal flow graph performance. 
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Introduction 

Nowadays, the dc-dc converters have mostly used in many applications such as a power supply in electronic 
systems. In the portable devices with the battery, the converter converts the input voltage to the proper output 
voltage to use it by the electronic load. The battery voltage can vary in a wide range depending on a charge level. 
At the low charge level, it may drop below the load voltage. Thus, to continue supplying the constant load voltage, 
the converter must be able to work in both step-up and step-down modes. In this paper, the Zeta converter is 
chosen to investigate. Zeta converter is a forth order dc-dc converter which is capable of bucking or boosting the 
input voltage. In order to provide an inside view of the converter's action, and obtain the needed information for 
controller design, modelling plays an important role [1]. 

Modelling dc-dc converters is an effective way to get open-loop frequency response and investigate the stability of 
the closed-loop systems. Modelling is the first stage for designing compensator or controller circuits. Previous 
works present different models for dc-dc converters [2]. Considering that the power electronic converters include 
nonlinear elements like switches and diodes, then modelling of them needs linearization. So, applying averaging 
and linearization techniques has a special importance. After linearization, the inside model of the system is 
extractable and therefore investigating the frequency response will be possible [3-5]. Modelling dc-dc converters by 
using small signal linearization and averaging techniques cause complexity in the obtained equations. Solving 
these equations for basic converters like; buck, boost and buck-boost is not a problem, but for high order 
converters, dealing with them will be more difficult [6]. In [7- 8], in order to derive a mathematical model of a dc- 
dc converter and study the transient states of it, a combination of Laplace and Z-transforms has been employed. A 
converter including switches is a nonlinear system, but it can be decomposed to two linear circuits; one for on-state 
and the other for off-state of the switch. These two linear circuits are illustrated by means of two signal flow 
graphs. Composition of two sub-graphs by using switching branches reaches to the graph of the whole converter. 
Switching branches are the only nonlinear parts of the converters. Thus modelling process is limited to the 
switching branches. In [9- 12], a signal flow graph method has been presented for modelling a dc-dc converter, 
then by using Mason's gain formula any desired transfer function of the system might be extracted and used to 
design any proper controller, but the presented graph method gets a sophisticated process to achieve the small 
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signal model. The proposed technique in this paper, gives greatly simplified mathematical and graphical 
representation of the system and is based on signal flow graph of the converter. This method has great advantage 
because of simplicity and being capable of giving any desired transfer function of the system to design a controller 
for a determined variable. Graphical modelling technique is used for surveying nonlinear and dynamic behavior of 
switching converters. This type of converter has different configurations and all have linear subsystems [13]. The 
aim of this paper is to derive the model of the Zeta converter by using the proposed signal flow graph based 
method. The open-loop system has some disadvantages. For example, the output cannot be compensated or 
controlled if there is a variation or disturbance at the input. For the case of Zeta converter, the changes in the input 
voltage or the load current will cause the converter's output voltage to deviate from the desired value. To design a 
proper controller or compensator a good model of the system is needed to investigate the stability by using 
different tools. One of these tools is frequency domain analysis [14-15]. In order to get frequency response of the 
converters, after modelling the root-locus plots and the frequency domain response in form of Bode magnitude and 
phase diagrams are achieved and analyzed. 

Many control strategies have been presented to obtain the desired output voltage in dc-dc converters. In past, the 
controller design was based on using small signal linearization. Linear PID and PI controllers of the converters are 
usually designed by using standard frequency response techniques based on small signal model of the converter 
by using linear control theories such as Ziegler-Nichol's method, root locus technique, hysteresis method. Bode 
plot and etc. The Ziegler-Nichols method is an experimental one that is widely used. 

The main contribution of this paper is to derive the model of the converter using the presented signal flow graph 
method and then designing a model- based controller for the Zeta converter. In this process, by applying the 
mentioned signal flow graph technique and Mason's gain formula, the input to output, control to output transfer 
functions are obtained. Then, the frequency domain responses in form of Bode magnitude and phase diagrams are 
achieved and the proper tuned PID controller is designed. At last, the simulation results based on the presented 
modeling and controlling methods are brought to show the good performance and behavior of the modeling and 
control methods. 



(a) The power circuit (b) Equivalent circuit in on-state. 



(c) Equivalent circuit in off-state. 
FIG. 1. ZETA CONVERTER. 


Review on DC-DC Zeta Converter 

The power circuit of the Zeta converter shown in Fig. 1(a) has one switch, one diode, two inductors and two 
capacitors. In this figure, y and v a are the source voltage and the output or load voltage, respectively. v a , v c and 
v Co are the voltages on capacitors C’ , C and C a , respectively. , and i L2 are inductors currents. The load is 
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supposed as a resistor, R . Zeta converter can perform to generate a positive voltage of the input voltage which can 
be less or higher than input. The converter acts in two operating modes, one of them is related to on-state and the 
other is the off-state of the switch. Figs. 1(b) and 1(c) illustrate the equivalent circuits of the on-state and off- state 
operating modes, respectively. 


Extraction of the State Equations 

Here, the open-loop performance characteristics of the converter are tested. The starting point is the extracting a 
mathematical and analytical model of the system. The model is a mathematical description of the behaviour of the 
real system that is enough for performing stability test. 

The state equations are extracted by using kirchhoff's voltage and current laws for the two cases of on and off 
states of the switch. These equations are calculated in loops and nodes with just one passive element, so all 
obtained equations are first order ones. One of the advantages of this method is in high order circuits, with 
increasing the order of the circuit and energy saving elements, the all extracted state equations will be of the first 
order and reduce the complexity of the equations and ease the analysis. The state equations are extracted for two 
operating modes. The first mode is when the switch is on and the diode is off, and the second mode refers to the 
case that the switch is in off-state and the diode conducts. The state equations for the circuit shown in Fig. 1(b) are 
as follows: 


dt 


T d ^L2 

L 2 —r- = Vi+v c -v 0 

dt 


dt L2 


v o - v Co 


C *0 


— IfJ 


dt 12 R 


Substituting Eq. (4) in Eq. (5), the following equation is achieved: 


c ^ = ; 

" ' ” R 


"° dt L1 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


The state equations for the case that the switch is off and the diode conducts (Fig. 1(c)) are as below: 

l f^ = _ v 

1 dt c 


■ = — V„ 


L d' L2 
dt 

dt u 
V o = V Co 

r dv co _ . 

0 dt 12 R 

Substituting Eq. (10) in Eq. (11), the following equation is achieved: 

0 dt L - R 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 
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Extraction of the Averaged State Equations 

The state equations at the on time interval of the switch is multiplied by d and at the off time interval is multiplied 
by 1 -d , then they are added to each other, d is considered as the switch conduction coefficient and 1 — d as the 
switch non-conduction coefficient, d can vary in 0 and 1 interval. 

In the following equations, the variables in form of x are average values which are different from the 
instantaneous variables in form of x . As mentioned above, Eq. (1) to Eq. (3) and Eq. (6) are multiplied by d , and 
Eq. (7) to Eq. (9) and Eq. (12) are multiplied by 1 — d . Then sum of the couples will be calculated and the result is 
the averaged state equations. Hence, the averaged equations are related to both on and off time intervals of the 
switch as follows: 


L \ —r~ = dv t - (1 - d)v c 
at 


T ^L2 J — 7 — — 

— — = d v t + d v c - v Q 


dt 


cA-\A 

0 dt L1 R 


,dVc_ —-r 
dt 


C — di L2 +( 1 d)i n 


(13) 

(14) 

(15) 

(16) 


Small Signal Linearization 

In averaged state equations, dc and ac variables are substituted and rewritten with X and x , respectively. 
Considering all variables in form of x = X + x , Eq. (13) to Eq. (16) are rewritten as below: 


L l s(i Ll +hi) = (D + d)(v i +V i )-(l-D-d)(v c + V c ) 


(17) 


L 2 s(i L2 +I L2 ) = (D + d)(v, + V,) + (D + d)(v c +V c )-v 0 -V 0 


CAv 0 + v 0 ) = i L 2 +i L 2 -^f-'f 

K K 


(18) 

(19) 


Cs(v c + V c ) = -(D + d )(i L2 + 1 L1 ) + (1 - D - d )(i Ll + I Ll ) 


( 20 ) 


The product of two dc variables is a dc one and the product of one dc and one ac variable is an ac variable, also the 
product of two ac variables is zero. Thus, Eq. (17) to Eq. (20) are obtained as below: 


L,siu=^,+(V,+V c )d-(\-D)v c 


( 21 ) 


L 2 sI l2 = DVj + (V. +V c )d + Dv c - v a 


C 0 sv 0 =i L2 —j- 

K 


( 22 ) 

(23) 


Csv c - Di L2 ( I L2 + I Li )d + (1 D)i L1 

Using dc equations, V c , /, and I L2 values versus V j and D are calculated as follows: 


1 -D 


(24) 


(25) 
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- — — DV 

V = DV t + DV C - — = 
1 -D 


1 - V sl- 

1 L2 


DV; 


R R{\-D ) 

dT l2 =(l-D)T u 


Il\ ~ , — I L2 ~ 


— — 2 
V;D 


1 -D " R(\-Dy 


(26) 

(27) 

(28) 
(29) 


Considering V c , I u and I L2 values and replacing them in Eq. (21) to Eq. (24) and simplifying them, the small signal 
ac equations are calculated as below: 


r D _ V t 1 -D 

si Ll = —V; H 1 


L \ ' U (1 - D) L x 


D 


V; 


D _ 1 


si T1 = — v, H ' _ d h v r 

L 2 L 2 (l-D) l 2 l 2 

1 ~ 1 _ 
sv o=-pr l Lz-—v 0 


D 


1 -D 


C RC(\-D) 2 C 


(30) 

(31) 

(32) 

(33) 


Proposed Signal Flow Graph 

The ac small signal equations given in Eq. (30) to Eq. (33) are used to plot the signal flow graph as shown in Fig. 2. 

In these equations, there are four variables; v a , i L] , i L2 and v c , and two inputs; d and v ; . The graph is depicted 
for all independent variables (sic, ic) by using nodes and paths between two nodes. In order to extract the transfer 
functions from signal flow graph, the Mason gain formula will be used. 



FIG. 2. SIGNAL FLOW GRAPH FOR ZETA CONVERTER. 


Extracting Transfer Functions Based on the Proposed Method 

For calculating the transfer functions, first the gain of distinct loops, non-touching loops and forward paths from 
each input to output are extracted as Table 1. 


Transfer functions v a /v i and v„ / tl are calculated by using Mason's gain formula youth in = £ Phk / A ■ In this 


formula, n is the number of forward paths from input to the output and A is calculated as follows: 
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A = 1 - X ( distinct loops gain) 

+ X ( double nontouching loops gain) 

- Z ( tripple nontouching loops gain)i — 

Omitting the path P n , A, =1 is obtained, and omitting the path P i2 , A 2 = 1 - Loop l is obtained. Replacing the 
obtained values, the transfer function v a / v. will be as shown in Table 1. Omitting the path P n , A, =1 is obtained, 
and omitting the path Pd2 , A 2 = 1 is obtained. Also, removing P d 3 , A 3 = 1 -Loop 1 is achieved. Replacing the 
obtained values, the transfer function v 0 / d will be obtained as given in Tablet. 
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TABLE 2. SIMULATION PARAMETERS 


A 

l 2 

C 

C, 

18.9 mH 

18.9 /uH 

14 nF 

11.2 nF 

c 0 

R a 

Vj 

d 

6.5 fiF 

ion 

12V 

0.57 


TABLE 3. INFORMATION OBTAINED FROM POLES AND ZEROS OF THE 
TRANSFER FUNCTION OF vjv, 



Poles 

Zeros 

Value 

— 1 . 1 5 x 1 0 3 
+2.45x1 Oh 

-6.54xl0 3 

+9.69xl0 4 ; 

+4.03 xlO 4 ! 

Damping 

0.047 

0.0673 

0 

Overshoot (%) 

86.3 

80.9 

100 

Frequency (rad/s) 

2.46 xlO 4 

9.71xl0 4 

4.03xl0 4 


X 10 


Root Locus 


-15000 


0 


5000 


-10000 -5000 

(a) ROOT LOCUS 

FIG 3. RESPONSE OF TRANSFER FUNCTION OF vjv, 



Simulation Results 

In simulation, the parameters of Zeta converter are considered as given in Table 2. 

Frequency Response Analysis for v a /v i Transfer Function 

The root-locus, magnitude and phase Bode diagrams for extracted transfer function of v 0 /v, are shown in Fig. 3. 
The information obtained from these diagrams is given in Table 3. 

Considering the root-locus shown in Fig. 3(a), it is understood that v D /v, has four complex poles in left hand side 
of the imaginary axis and two complex zeros on imaginary axis. Considering Bode diagrams (Fig. 3(b)), the 
magnitude curve of the function begins with a constant value of 2.64 dB (DC gain) and at the frequency of 24 ,6kHz 

duo to the poles -1.15xl0 3 ±2.45xl0 4 ; , an undershoot occurs in magnitude curve and its slope becomes 
-40 dB / dec and the phase curve begins from 0° and a phase reduction of -180° takes place. At the frequency of 
40.3 kHz duo to the zeros +4.03 x 10 4 / on imaginary axis, an overshoot occurs in magnitude curve and +40 dB I dec 
is added to its slope. So, the slope becomes zero and +180° is added to phase curve duo to the mentioned zeros. 
Finally at frequency 97.1 kHz because of the poles -6.54 x 10 3 ± 9.69 x 1 0 4 / , an undershoot occurs in magnitude curve 
and its slope becomes -40 dB / dec and at the phase curve a phase reduction of -180° takes place and the phase 
becomes -180° . 

Considering the plotted responses for all values of the gain, the transfer function v Q /v, is stable and gain margin 
becomes co. Considering Bode diagram, when the magnitude is 0 dB , the phase equals to 171°, so the phase 
margin which equals to the difference of 180° and this phase will be 9° . 

2) Frequency Response Analysis for v a /d Transfer Function 

The root-locus and Bode diagrams for the extracted transfer function of v a jd are shown in Fig.4. The information 
obtained from these diagrams is given in Table 4. 

Considering the root-locus shown in Fig. 4(a), it is understood that v 0 jd has four complex poles in left hand side 
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of the imaginary axis and two complex zeros right hand plane. From Bode diagrams (Fig. 4(b)), gain margin of the 
system is -37.5 dB and the phase margin is 2.52° . So, the transfer function is unstable and needs a compensator in 
normal conditions. 


From root-locus diagrams, the mentioned system has two complex zeros in right hand side of the imaginary axis. 
In order to have a more stable system the proper controller should be designed for the mentioned transfer function. 
In next section designing the proposed model- based PID controller will be explained. 


x 10 


Root Locus 


4 

2 

0 - 
-2 

-4 - 



-1.5 -1 -0.5 0 

(a) ROOT LOCUS 

FIG 4. RESPONSE OF TRANSFER FUNCTION OF vjd 



TABLE 4. INFORMATION OBTAINED FORM POLES AND ZEROS OF THE TRANSFER FUNCTION OF vjd 



Poles 

Zeros 

Value 

-l.I5xI0 3 
±2.45 xIO 4 ; 

-6.54 xIO 3 
+9.69xl0 4 i 

2.7 xIO 3 
±4.03 x10 4 i 

Damping 

0.047 

0.0673 

-0.0669 

Overshoot (%) 

86.3 

80.9 

123 

Frequency (rad/s) 

2.46 xIO 4 

9.71xl0 4 

4.03 xIO 4 


Step Response 




FIG. 6. THE WHOLE SYSTEM CONFIGURATION WITH THE 
PID CONTROLLER 


The Proposed Model-Based Controller 

Fig. 5 shows the step response of the system. This shows that the open-loop system does not reject disturbances on 
the input voltage and cannot regulate the output voltage. 

At the following, a model-based PID controller by using Ziegler-Nichol's tuning method will be explained. 

The whole configuration of the controller beside the controlled system is plotted in Fig. 6. In order to control the 
converter based on the previously extracted model, the PID controller is added to the controlled transfer function 
of the converter. 

The PID controllers are widely used because they have only three parameters that has to be tuned for the process 
control. The input of a PID controller is an error signal which is the difference between the measured value and a 
desired reference signal for a process variable. The controller minimizes the error by tuning three constant 
parameters namely the proportional term (K p ) , the integral term (K t ) and the differential term ( K D ) . Because of 
its simple configuration, the PID controller is the most commonly used control method in industry. The important 
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aims of the PID controllers are: 

> eliminating steady-state error of the step response (due to the integral action) 

> reducing the peak overshoot (providing damping due to derivative action) 

There are many methods for tuning the PID parameters or gains. The most famous method for achieving the 
mentioned goal is Ziegler-Nichols tuning method. Table 6 shows the parameters obtained for the controller. 

Fig. 7 shows the step response of the system after applying the designed PID controller. 


Fig. 8 shows the Bode diagrams of the system before and after applying the designed PID controller. Fig. 9 shows 
the Nyquist diagram of the closed loop system. The results prove the stability of the obtained system with the 
proposed controller. Stability of the obtained controlled system is well illustraited in Fig. 9. 


From the simulation results, it is obvious that the disturbance is rejected and this controller behaves very desirably 
in time-domain. Checking the bode diagrams depicts that the system has become a stable one. Also, calculations 
bring out the gain margin as 16.2 dB and a phase margin of 89.7 deg which are favorite frequency-domain 
response characteristics. Then, applying the proposed controller, the converter system becomes stable and 
regulated. 


TABLE 6. PID PARAMETERS OBTAINED 


K P 

0 

K, 

8.0306 

k d 

0 

Rise time of dosed loop response 

0.00422 

Settling time of closed loop response 

0.00748 

Overshoot of dosed loop response 

0 

Gain Margin of dosed loop response 

16.2 

Phase Margin of dosed loop response 

89.7 

Closed loop stable? 

Yes 


Step Response 



FIG. 7. STEP RESPONSE OF THE CLOSED LOOP CONVERTER 
SYSTEM 


Bode Diagram 



FIG. 8. BODE DIAGRAMS OF THE OPEN LOOP AND 
CLOSED LOOP CONVERTER SYSTEM 
Nyquist Diagram 



FIG. 9. NYQUIST DIAGRAM OF THE CLOSED LOOP 
CONVERTER SYSTEM 


Conclusions 

In this paper, modelling and analyzing of the Zeta converter operating in continuous conduction mode were 
presented. Application of the proposed based on signal flow graph method and Mason's gain formula for 
extracting associated transfer functions and modelling a dc-dc Zeta converter was investigated. In this process, for 
modelling the average and small signal linearization techniques were employed. Extracting the final transfer 
functions from initial state space equations were described step by step. Application of signal flow graph method 
and Mason's rule helps to solve the high order complex equations simply. The presented method deals with the 
equations with more state variables as well. In this paper, the transfer functions from input to output and control to 
output were obtained. Then, by using the proposed method the poles, zeros, root-locus plot and Bode diagrams 
were extracted and the achieved frequency response was analyzed. Finally, the efficiency of the proposed tuned 
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PID controller as a model based technique to design a proper control system is shown. The proposed control 

system employs the PID controller with a special adjustment for proper setting of the control responses. 
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